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METHOD AND APPARATUS FOR EXPANDING A COLOR GAMUT 



Field of the Invention 



This invention relates to color printing, and more particularly, to the use of 
colorants and gamut scaling in color printing. 



5 



Background of the Invention 



Combinations of the colorants cyan, magenta, and yellow (CMY) are capable of 
producing images containing colors in the visible spectrum, so many color imaging 
systems are designed around CMY colorants. To improve image quality and to reduce 
the costs of printing, black is sometimes added to a CMY color gamut to form a CMYK 

10 (cyan, magenta, yellow, and black) color gamut. The visual details of an image rendered 
from a CMYK color gamut are sharper than the visual details of the same image rendered 
from a CMY color gamut. Similarly, the gray balance of an image rendered from a 
CMYK color gamut is better than the gray balance of an image rendered from a CMY 
color gamut. So, many modern color imaging systems utilize the colorants cyan, 

15 magenta, yellow, and black (CMYK) in forming color images. 

Most conversion methods from CMY to CMYK simultaneously add black 
colorant and reduce the amounts of cyan, magenta, and yellow. One popular solution is 
based on a combination of Under Color Removal (UCR) and Black Generation (BG) 
functions: 



25 where (C, M, Y) are input values and (C 1? M„ Y 1? K, ) are the output values. The black 



amount P. The CMY values are reduced by an amount dependent on PK. However, such 



20 



PK = min(C, M, Y) 
C, = C - UCR(PK) 
M, = M - UCR(PK) 
Y, = Y - UCR(PK) 
K,=BG(PK) 



(1) 



generation function computes the black colorant amount based on the process black 
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a simple transformation does not preserve colorimetric relations and does not allow use of 
the whole printer gamut. 

Unfortunately, even though each of these processes show improved gray balance 
in images rendered from a CMYK color space when compared with images rendered 
5 from a CMY color space, each of these processes also has some disadvantages. One 
disadvantage is that the UCR and GCR processes are not based on human visual 
perception and therefore fail to faithfully reproduce the original image. A second 
disadvantage of most CMYK systems is that they fail to utilize the entire gamut of the 
CMYK color space. Failing to utilize the entire gamut results in a suboptimal rendering 
10 of a CMY image in a CMYK color space. 



The above mentioned problems with generating a CMYK color gamut from a 
CMY color gamut are addressed by the present invention, which will be understood by 

15 reading and studying the following specification. 

In one embodiment, a gamut in a first color space is converted to a gamut in a 
second color space. The gamut in the second color space is converted to a gamut in a 
third color space. A lightness component of a gamut value in the third color space is 
rescaled to form a modified gamut. Finally, the modified gamut is converted to a second 

20 gamut in the second color space. In an alternate embodiment, rescaling a lightness 

component of a gamut value in the third color space to form a modified gamut comprises 
modifying the gamut in the third color space by changing a lightness component of the 
gamut value in the third color space such that an upper surface of the gamut in the first 
color space is preserved and a lower surface of the gamut in the first color space is 

25 mapped to a bottom surface of the gamut of the second color space to form an expanded 
gamut in the third color space. 

In another embodiment, a gamut in a CMY color space having an upper surface 
and a lower surface is converted to a gamut in a CMYK color space having a bottom 
surface. The gamut in the CMYK color space is converted to a gamut in a CIELAB color 



For these and other reasons there is a need for the present invention. 



Summary of the In ve ntion 
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space and the gamut in the CIELAB color space has a lightness component. The gamut in 
the CIELAB color space is modified by changing the lightness component such that the 
upper surface of the gamut in the CMY color space is preserved and the lower surface in 
the CMY color space is mapped to the bottom surface of the gamut of the CMYK color 
5 space to form a gamut in an expanded CIELAB color space. Finally, the gamut in the 
expanded CIELAB color space is transformed to form a second gamut in the CMYK 
color space. In an alternate embodiment, a CIELAB gamut having an L* component is 
modified to form a modified CIELAB gamut, and the modified CIELAB gamut is 
converted to a color control space gamut subject to K constraints and Total Area 
10 Coverage (TAC) constraints. In another alternate embodiment, the CIELAB gamut is 

modified to form a modified CIELAB gamut and the method comprises linearly rescaling 
the L* component of the CIELAB gamut to form the modified CIELAB gamut. 
=£ In another embodiment, a method comprises obtaining a CIELAB space gamut 

'jph from a CMYK space gamut obtained from a CMY space gamut, expanding the CIELAB 

Vsb" 

]}f s 15 space gamut by changing the lightness component of the CIELAB space gamut, and 
Ijs transforming the CIELAB space gamut into a second CMYK space gamut. In an 

'i** alternate embodiment, changing the lightness component of the CIELAB space gamut 

comprises linearly rescaling the lightness component of the CIELAB space gamut. In 
another alternate embodiment, linearly rescaling the lightness component in the CIELAB 
20 space gamut comprises computing a rescaling factor that is a function of an L min , an L max , 
and an L mincmy . 

In another embodiment, an apparatus for transforming a CMY color space gamut 
into a CMYK color space gamut comprises a processing unit, a memory unit coupled to 
the processing unit, and a software means operative on the processing unit. The software 
25 means is operative on the processing unit for obtaining an CIELAB space gamut from a 
first CMY space gamut. The CIELAB space gamut has a plurality of values and each of 
the plurality of values has a lightness component. After obtaining a CIELAB space 
gamut from a first CMY space gamut, the software means is operative for expanding the 
CIELAB space gamut by changing the lightness component of at least one of the plurality 
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of values of the CIELAB space gamut and transforming the CIELAB space gamut into 
the CMYK space gamut. 

In another embodiment, a computer-readable medium has computer-executable 
instructions for performing a method. The method comprises first obtaining a CIELAB 
space gamut having a lightness component. The CIELAB space gamut is obtained from a 
CMYK space gamut which is obtained from a CMY space gamut. Second, the method 
comprises expanding the CIELAB space gamut by linearly rescaling the lightness 
component in the CIELAB space gamut using an L min , an L max , and an L mincmy Finally, the 
method comprises transforming the CIELAB space gamut into a second CMYK space 
gamut. 



Figure 1 is a flow diagram of example embodiments of a method for black 
generation using lightness scaling. 

Figure 2 is a diagram and a visualization of the lightness rescaling process in the 
CIELAB space. 

Figure 3 is a projection for a color inkjet printer gamut surface. 

Figure 4 is a diagram of alternate example embodiments of a method for 
increasing the use of a color gamut. 

Figure 5 is a high level diagram of one embodiment of the inverse modeling 
process shown in Figure 4. 

Figure 6 is a block diagram of example embodiments of an apparatus for 
transforming a CMY color space gamut into a CMYK color space gamut. 

Figure 7 is a diagram of one embodiment of a computer readable media including 
computer executable instructions. 



The present invention provides a method for relating a three-dimensional CMY 
coordinate system to CMYK colorant amounts in a way that increases the utilization of a 
printer gamut volume. A printer gamut, in one embodiment, is an arbitrary set of colors. 



Brief Description of the Drawings 



Detailed Description of the Invention 
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In general, adding a black colorant to any combination of CMY colorants does not result 
in a color with more chroma. Therefore, a CMYK gamut can be obtained from a CMY 
gamut by expanding the CMY gamut towards lower lightness values. Use of black 
colorant on the gray axis is enforced by modifying the initial distribution of CMY points 
5 through an approximate black generation transform. Lightness values of a resulting set of 
points in CIELAB space are scaled to fill the four-color gamut volume. The CMYK 
values corresponding to the modified CIELAB colors are found by inverting a printer 
model. This last process determines a specific black use rate which can depend on the 
region of the color space. 

10 Figure 1 is a flow diagram of one embodiment of method 100 for black generation 

using lightness scaling. Black generation method 100 converts a CMY color space to a 
CMYK color space and includes black approximation process 103 for converting the 
CMY color space to a CMYKj color space. Black approximation process 103 is followed 
by forward model process 106 for converting the CMYKj color space to a CIELAB color 

15 space. The CIELAB color space is converted by lightness scaling process 109 to a Lab! 
color space. Finally, inverse model process 1 12 converts the Lab! color space to a 
CMYK color space. 

In black approximation process 103, an approximate conversion from a CMY 
color space to a CMYK, color space is performed. Colors corresponding to output values 

20 generated by this conversion are within the original CMY gamut chroma boundaries. In 
addition, color values corresponding to the upper surface of the printer gamut are left 
unmodified, i.e.: 



In one embodiment of black approximation process 103, a K=0 component is 
25 added. In an alternate embodiment, to impose a condition that balanced CMY values be 
converted to black: 



Ci = C, Mi = M, Y, = Y, Ki = 0, if min(C, M, Y) = 0 and K = 0. 



C!=M! = Y, = 0, Ki=C, ifC = M = Y. 
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All the above conditions are satisfied by a black generation method known as 
achromatic reproduction, which is obtained from Equation (1) when BG(x) = x and 
XJCR(x)=x. 

K, = min(C, M, Y) 
5 C, = C - K, 

M, = M - K, 
Y,=Y-K,. 

The final a* and b* coordinates of colors on ramps from black to primary and 
secondary colors are also implicitly defined in black approximation 103. By using 

10 nonlinear modifications of the above formula, desired colorant behavior can be obtained. 

In one embodiment of forward model process 106, CIELAB values corresponding 
to the CMYKj values are obtained by printing and measuring a set of patches on a printer. 
In an alternate embodiment, CIELAB values are obtained from the CMYK, values by 
processing the CMYK, values using a printer model. A printer model is flexible and 

15 avoids the time and expense required to print and remeasure the patches when the black 
function is changed. A number of different techniques can be used for implementing a 
printer model. For example, a printer mode can be implemented as Neugebauer 
equations, lookup tables with interpolation, neural networks, fuzzy logic, or a polynomial 
regression. 

20 In lightness rescaling process 109, the CIELAB point distribution is modified by 

changing the lightness components in such a way that the upper surfaces of the CMY 
gamut are preserved and the lower surfaces are mapped to the bottom surface of the full 
CMYK gamut. For each point obtained in forward model process 106, three 
corresponding points with identical a* and b* coordinates are identified. One of the points 

25 is located on the top surface of the CMYK gamut, which is identical to the top surface of 
the CMY gamut. The lightness of this point is denoted by L* max . The other two points are 
located on the bottom surfaces of the CMY and CMYK gamuts and have lightness 
L*mincmy an d L* min , respectively. If black is not the darkest point in the gamut, to avoid 
non-monotonic behavior in some regions of the transform, the values of L* min can be 
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clipped to not less than L* of black. In one embodiment, the modified L* values are found 
by expanding the lightness range by linear rescaling: 



Lx = 



^min + ~ ^mincmy) C^max ~ ^min) ^ C^max ~ ^mincmy) " ^max ~ ^mirtcmy * ^ 

I * otherwise. 



A diagram and visualization of the lightness rescaling process is shown in Figure 2. 
In one embodiment of the present invention, a geometric method is used to 
5 determine the L* min , L* max and L* mincmy values. In an alternate embodiment, a constrained 
optimization method is used to determine the L* mjn , L* max and L* mincmy values. Neither of 
these embodiments is based on assumptions about the way CMY sub-gamuts for different 
values of K overlap and work in the presence of natural gamut boundaries. 

The geometric method uses alpha shapes to determine the description of the 

10 gamut boundaries in the form of surfaces composed of triangles. The intersections of the 
lines of constant (a*, b*) with the top surface of the CMYK gamut and the bottom surfaces 
of CMYK and CMY gamuts are found. A search for triangles which are intersected by 
this line is performed after the L* coordinates are ignored by projecting all the triangles 
onto the aV plane. Figure 3 shows such a projection for a color inkjet printer gamut 

15 surface. The solid and dashed lines correspond to the top and bottom gamut surfaces, 
respectively. Planar barycentric coordinates of a point (a*, b*) with respect to a triangle 
P 1 P 2 P 3 where P ( = (L * , a* , b* ) are given by: 

u = [(a 2 - a ) (b* 3 - b') - (a* 3 - a*) (b* 2 - b*)] I A 
v = [(a 3 - a ) (b* t - b*) - (a, - a ) (b* 3 - b*)] I A 
20 w = [(a\ - a ) (b* 2 - b*) - (a 2 - a*) (b\ - b*)] I A 

where A is the area of the triangle under consideration: 

A = (a* 2 - a*,) (b* 3 - b*,) - (a* 3 - a*,) (b* 2 - b\). 
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The line of constant (a*, b*) intersects a triangle if each of the coordinates (w, v, w) 
is greater than or equal to zero. Assuming that the gamut shape is sufficiently regular, for 
each (a*, b*) within the gamut chroma limits there will be two triangles like that: one on 
the top and one of the bottom surface of the gamut. The L* value for each of the 
5 intersection points is then computed from the weighted sum: 

L* = uL^ + vL 2 + wL v 

Most printing processes impose a limit on the total Area Coverage (TAC), defined 
as the sum of the fractional areas covered by the halftone dots in the four separations. 
Exceeding the TAC limit is avoided by using only the points with colorant amounts less 

10 than TAC for computation of gamut surfaces. 

The constrained optimization method is adapted for the purpose of finding the 
maximum and minimum L* values on the gamut surface by searching for the point in the 
gamut which is the closest to the point outside the gamut which has the same a* and b* 
coordinates. The color space error function is constructed by weighting the a* and b* (or 

15 chroma and hue) much more strongly than lightness, for example, by setting S c = S H = 1 
and S L = 0.001 . This produces the type of gamut mapping which transforms out-of- 
gamut points approximately vertically to points on the gamut surface. The range of 
lightness values is limited by the total area coverage (TAC) constraint incorporated in the 
optimization process. 

20 For inverse modeLprocess 112, the expanded CIELAB values are inverted. The 

inversion process is subject to specific constraints. One constraint is based on the desired 
relationship between K and CMY. A second constraint is based on the Total Area 
Coverage (TAC) limit. The TAC constraint allows the black use rate to vary with the 
region of the color space to which the inversion is being applied. 

25 After the CMYK values are obtained by inverse model process 112, the CMYK 

values can be smoothed. One useful smoothing kernel is based on a Laplace equation, 
which in the one-dimensional case (for the edges) is: 
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Uj = - (Uj_ x + u J+l ) for ; = 1,. . . , J - 1 , 



where J is the number of grid points. In one embodiment, the smoothing process includes 
three operations. First, CMYK values for edges of the input CMY cube are smoothed, 
while keeping the corner values fixed. Second, the cube faces of the CMY cube are 
smoothed, while keeping the edges fixed. Finally, the whole cube volume is smoothed. 
5 Alternatively, smoothing can be performed by fitting a global model to the set of CMYK 
values and then using the fitted global model to obtain new values. 

Figure 2 is a diagram and a visualization of the lightness rescaling process in the 
CIELAB space. Dashed line 203 shows the bottom of the CMY gamut. Line segments 
206 and 209 show the intermediate outline of the gamut after approximate black 

10 generation. Line segments 212 and 215 show an outline of the full CMYK gamut. 
Vertical dashed line 218 has constant a* and b* coordinates. Point P t is not modified 
because it is on the upper surface of the gamut. Point P 2 is mapped from the bottom of 
the modified CMY gamut to a point P 3 on the bottom of the CMYK gamut. Point P 4 is 
inside the CMY gamut and is transformed to a point P 5 which has proportionally rescaled 

15 lightness. 

Figure 4 is a diagram of alternate example embodiments of method 400 for 
increasing the use of a color gamut. The diagram includes first color space gamut 403, . 
second color space gamut 406, third color space gamut 409, black generation 412, gamut 
distortion 415, forward mapping 416, linear rescaling 418, inverse modeling 419, 

20 modified gamut 421, and second gamut 424. The diagram also includes conversion paths 
427, 430, 433, 436, and 439. 

In one embodiment, first color space gamut 403 is a CMY color space gamut, 
second color space gamut 406 is a CMYK color space gamut, and third color space gamut 
409 is a CIELAB gamut. In one embodiment of method 400, conversion path 427 

25 employs black generation 412 in converting first color space gamut 403 to second color 
space gamut 406. In an alternate embodiment, conversion path 430 employs geometric 
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distortion 415 in the process of converting first color space gamut 403 to second color 
space gamut 406. 

After converting first color space gamut 403 to second color space gamut 406, 
conversion path 433 employs forward mapping 416 to convert second color space gamut 
5 406 to third color space gamut 409. In one embodiment, forward mapping 416 includes 
printing color patches from second color space 406 and performing a spectral analysis of 
the patches. The spectral analysis provides values in third color space gamut 409 into 
which the values of second color space gamut 406 may be mapped. In an alternate 
embodiment, forward mapping 416 includes developing and employing a model to map 
10 second color space 406 to third color space 409. 

After converting second color space gamut 406 to third color space gamut 409, 

2 conversion path 436 employs linear rescaling 418 to convert third color space gamut 409 

y 

P to modified gamut 42 1 . 

ig Finally, after converting third color space gamut 409 to modified gamut 421, 

15 conversion path 439 employs inverse modeling 419 to convert modified gamut 421 to 
y\ second gamut 424. 

^ Figure 5 is a high level diagram of one embodiment of inverse modeling 419 

shown in Figure 4. A forward model f(x) 503 for mapping a control space 506, such as a 
CMYK color space, to a target space 509, such as a CIELAB color space is developed. 
20 An inverse model g(x) 5 12 for mapping the target space to the control space is derived 

from the forward model and various constraints. This method permits the specification of 
constraints on the desired relationship between K and CMY and on the Total Area 
Coverage (TAG) limit. This conversion process also produces second gamut 424 from 
modified gamut 421 as shown in Figure 4. 
25 Figure 6 is a block diagram of example embodiments of apparatus 600 operable 

for transforming a CMY space gamut into a CMYK space gamut. Apparatus 600 
includes processing unit 603, software means 606, storage unit 609, and image forming 
unit 612. Processing unit 603 is coupled to storage unit 609 and image forming unit 612. 
Software means 606 is operative on processing unit 603. Color image 615 is an input to 
3 0 processing unit 603 . 



11 
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Processing unit 603 is not limited to a particular type of information processor. 
Processing units suitable for use in connection with the present invention include 
microprocessors, mainframe processors, mid-range processors, reduced instruction set 
processors, complex instruction set processor, parallel processors, and massively parallel 
5 processors. 

Storage unit 609 is not limited to a particular type of storage unit. In one 
embodiment, storage unit 609 is a solid state memory. Alternatively, storage unit 609 is a 
magnetic memory. In still another alternate embodiment, storage unit 609 is an optical 
memory. 

10 Image forming unit 612 is not limited to a particular type of image forming unit. 

In one embodiment, image forming unit 612 is a color image forming unit, such as a 
color ink jet printer, a color laser printer, or an offset printing system. Alternatively, 
C image forming unit 612 is a display device, such as a CRT, a plasma display, or a liquid 

g crystal display. In still another alternate embodiment, image forming unit 612 is an 

j? 1 5 image capture device, such as a color scanner. 

T- Software means 606, in one embodiment, performs the operations of method 100 

for black generation using lightness scaling shown in the flow diagram of Figure 1. 

;* Figure 7 is an illustration of one embodiment of computer-readable medium 700 

including computer executable instructions 703. Computer readable medium 700 

.Jj 

20 includes any media capable of storing information, such as magnetic media, optical 

media, or electronic semiconductor media. Computer executable instructions 703 include 
instructions for expanding a color gamut, such as could be coded based on the method 
shown in Figure 1 . 

A method and apparatus for expanding a color gamut has been described. The 
25 method includes a series of conversions and transformations between color space gamuts 
which result in expanding and transforming an original color space gamut in such a way 
that a control space gamut, such as a CMYK space gamut, is more completely utilized. 

Although specific embodiments have been illustrated and described herein, it will 
be appreciated by those of ordinary skill in the art that any arrangement which is 
30 calculated to achieve the same purpose may be substituted for the specific embodiment 
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shown. This application is intended to cover any adaptations or variations of the present 
invention. Therefore, it is intended that this invention be limited only by the claims and 
the equivalents thereof. 
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